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The Pyragas method for controlling chaos is investigated 
in detail from the experimental as well as theoretical point 
of view. We show by an analytical stability analysis that 
the revolution around an unstable periodic orbit governs the 
success of the control scheme. Our predictions concerning 
the transient behaviour of the control signal are confirmed by 
numerical simulations and an electronic circuit experiment. 
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I. INTRODUCTION 

The problem of controlling unstable motion is a clas- 
sical subject in engineering science. The revived inter- 
est of physicists in this subject, however, started with 
the observation that a large number of unstable periodic 
orbits embedded in chaotic attractors can be stabilized 
by weak external forces Q. Since that time a real in- 
dustry on chaos control has developed [^J. Two main 
methods for controlling unstable motions have been es- 
tablished meanwhile. The first one, developed by Ott, 
Grebogi, and Yorke 0, is based on the invariant mani- 
fold structure of unstable orbits. It is theoretically well 
understood, but difficult to apply to fast experimental 
systems. The second one, proposed by Pyragas ||, uses 
time-delayed controlling forces. In contrast to the former 
one it can easily be applied to real experimental situa- 
tions, but so far the control mechanism has been poorly 
understood from a theoretical point of view. By perform- 
ing an analytical linear stability analysis we demonstrate 
which class of orbits is accessible to time-delayed feed- 
back control methods. In addition, we obtain explicit 
expressions for important quantities like the critical and 
optimal control amplitude or the dependence of the tran- 
sient behaviour on the control parameters. 

II. THEORETICAL APPROACH 

We consider a dynamical system which is described by 
a general set of differential equations. It may contain a 
periodic explicit time dependence. 
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x=ftp(t),t) (1) 

We are interested in the stabilization of an unstable pe- 
riodic orbit£(£) =£(t + t). t is an integer multiple of the 
period of the driving force for non-autonomous systems. 
We remind the reader that the linear stability analysis of 
such an orbit according tox{t) =£(t) +exp [(A + ioj)t)u(t) 
leads to a Floquet problem, where the exponent and the 
periodic eigenfunctionu(i) =u(£ + r) are determined by 

[X + ico^(t)+d = If^t),t]i(t) , (2) 

with Df denoting the Jacobian matrix off. Since our sub- 
sequent analysis applies separately to each Floquet expo- 
nent we refrain from numbering the different branches. 
The real and the imaginary parts of the Floquet expo- 
nents govern the instability and the revolution of the tra- 
jectory around the unstable periodic orbit (cf. fig.0). 




FIG. 1. Diagrammatic view of a trajectory in the vicinity 
of an unstable periodic orbit. 

In order to achieve control of the unstable periodic or- 
bit the system ([!]) is, following the idea of j|, subjected 
to a time-delayed force Kip(t) —x(t — r)]. The most gen- 
eral situation is given by 

z=F^(t),K^t)-z(t-T)],t) , (3) 

where the right-hand side obeys the constraintF|:, 0, t) = 
f^,t), and the amplitude K of the controlling force is 
introduced for convenience. As long as the delay time 
coincides with the period of the unstable periodic orbit 
the controlled system admits the same solutions (i) =£(t). 
Linear stability analysis according toz(t) + &(t) 

yields 

& = DFZ(t),Q,t)&(t) 

+ Dft(t),Q,t)K[&(t)-&(t-T)] , (4) 
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where Df denotes the Jacobian matrix with respect to 
the i th (vector type) argument. In the case of conven- 
tional Pyragas control, where only one system variable is 
assumed to be accessible, the matrix Df contains only 
one non-vanishing element on the diagonal. But we keep 
our approach as general as possible. The (infinite di- 
mensional generalization of) Floquct theory tells us 
that the deviations obey &(t) = exp [(A + iil)t]t}(t) and 
v(t) =v(t + r), so that eq.(Q) reduces to 

[A + iO]u(i) +v 
= A[K(l-exp[-AT-iClT]),t 1 p(t) . (5) 

Here the abbreviation 

A[ K ,t) :=#£(*),*) +I«P£(t),0,t)K (6) 

has been used. From eq.(||) it is obvious that A + ifl can 
be expressed in terms of the Floquet exponents of the 
matrix (|^). If we denote the latter for convenience by 
r [k] then eq. (||) implies the relation 

A + iQ = T[K(1 - exp [-At - iftr])] . (7) 

This expression, which in fact is not entirely new but 
has been evaluated only numerically for specific examples 
(cf . ||[| ) , determines the exponents of the controlled or- 
bit in dependence of the control amplitude K. Although 
it is in general a difficult task to obtain a closed analyti- 
cal expression for the quantity T, we know by definition 
that the boundary condition (cf. cqs.(^|) and (||)) 

T[0] = X + iLU (8) 

holds, and that T is an analytical function as long as 
the Floquet exponents are non-degenerate. These prop- 
erties are sufficient to conclude that only orbits with a 
finite frequency Q ^ can become stable. On increas- 
ing the control amplitude K the real part of the Floquet 
exponent A has to change its sign from positive to neg- 
ative values in order to achieve stabilization. But if the 
frequency Q of the controlled orbit remains zero, the in- 
fluence of the controlling force K[l — exp(— At)] vanishes 
if the orbit tends to become stable, so that the solutions 
A of eq.(0) never can change their sign. The reader might 
object that the condiditon SI = is atypical and does not 
occur generically. But we remind of the fact that (nondc- 
generate) real Floquet multipliers are stable with respect 
to perturbations (cf. 0), so that both cases $1 = and 
Q ^ occur in a sense with equal proability. However 
one should keep in mind that the neccessary condition 
fl ^ for stabilization has to be fulfilled for each Floquct 
branch separately. Hence stabilization may be unlikely if 
the unstable manifold is high dimensional. 

We summarize that by the Pyragas method only orbits 
with a finite torsion can be stabilized, since for stabiliza- 
tion the influence of the controlling force has to be finite 



(cf. fig.[y). This property has been observed recently even 
in high-dimensional dynamical systems by analysing the 
transient behaviour of the control signal ||], but no ex- 
planation has been proposed. From this point of view 
the control methods by Ott, Grebogi, and Yorke on the 
one hand and by Pyragas on the other hand are com- 
plementary, since the former is not in principle but in 
most practical applications restricted to the case of only 
one unstable eigendirection, whereas the latter requires 
at least a two-dimensional unstable manifold. 

For further quantitative investigations some informa- 
tion about r[/c] is required. There are a few cases where 
by inspection the K-dependence can be read off from 
eq.®: 

AF£(t),0,<) = 1 

=> T[k] = A + ioj + k (9) 

DFt(t),o,t) = m(t),t) 

=>T[k} = {X + iuj)(l + k) (10) 

But we do not intend to confine our analysis to these 
special situations (cf. ||). Instead we suppose that the 
controlling force is small enough in order to neglect higher 
order terms in the expansion of eq. (0) . 

A + iti = A + iu 
+ ( X ' + tx")K (1 - exp [-At - iOr]) + 0{K 2 ) (11) 

Here use has been made of relation (||), and the abbre- 
viation x' + ix" '■— dT/dn\ K= o contains the details of the 
coupling mechanism of the controlling force. It is worth 
to mention that relation (|Il|) is exactly valid for the cases 
described by eqs.(^), (|io|). 

Relation ([ll]) determines the stability of the controlled 
orbit in terms of the control amplitude K, the Floquct 
exponent of the uncontrolled orbit X + iuj, and the precise 
mechanism of the coupling x\ x" ■ 

A = A + K X ' (1 - cxp(-At) cos(fh-)) 

- K X " exp(-AT) sin(^T) (12) 

Q = lj + Kx' exp(-AT) sin(fh-) 

+ K X " (1 - exp(-AT) cos(^t)) (13) 

For the evaluation we confine the subsequent discussion 
to an uncontrolled unstable periodic orbit which just flips 
its neighbourhood within one period, that means to an 
orbit of frequency u = tt/t. Such a situation appears par- 
ticularly in a neighbourhood of a period doubling bifur- 
cation. Since the corresponding Floquet exponent X + iuj 
is located at the "boundary of the Brillouin zone", that 
means the corresponding multiplier is an isolated nega- 
tive real number, T[k] —ico is by definition a real function 
at k — (cf. eq.(||)), and x" vanishes. For this reason 
fi = 7t/t is a solution of eq.(p^) which in that sense is 
an optimal value for the frequency since the effect of the 
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control term in eq.(|l2|) becomes maximal. Then eq. (|Ts 
simplifies to 



A = A + Kx' (1 + exp(-Ar)) , ft = tt/t 



(14) 



Stabilization is achieved if A changes its sign, that means 
at a critical control amplitude 



K c = -A/(2 X ') 



(15) 



K c is mainly determined by the instability of the uncon- 
trolled orbit. On further increase of the control ampli- 
tude the frequency may start to deviate from its optimal 
value. Formally this deviation results from a pitchfork 
bifurcation in eq.([l3|) which occurs at K 



opt 



1 = -if op ix'Texp(-A opt T) 



with A opt being determined by eq.(|lj). Beyond K, 



(16) 
the 



Op I 

frequency Q deviates from its optimal value tt/t so that 
the eigenvalue A of the controlled orbit starts to increase 
again with the control amplitude K . In this sense K opt 
is the optimal value since the stability of the controlled 
orbit is maximal. 



III. SIMULATIONS AND EXPERIMENT 

In order to illustrate our theoretical considerations and 
to demonstrate that the features predicted are accessible 
from observed data, we have performed numerical simu- 
lations of the driven Toda oscillator 

z\ = z 2 

Z2 = — \iZ2 — a (cxp(zi) — 1) + Asm2iTt 

- K [z 2 (t) - z 2 (t - 1)} . (17) 

Time is measured in periods of the driving force. At 
/i = 0.8, a = 25, A = 105, and K = the system pos- 
sesses a chaotic attractor. A period-one orbit, which has 
become unstable in a period doubling bifurcation, can be 
stabilized at finite control K > K c m 2.1. From the expo- 
nential decay of the control signal z 2 (t) — z 2 (t — 1), using 
a standard least square method, we determine the real 
and imaginary parts of the Floquet exponent A + id for 
several values of the control amplitude. Our findings are 
summarized in fig.^[ We clearly observe the predicted 
dependence on the control amplitude. Using the value 
K opt ~ 2.4 we determine the two parameters A and x' 
by means of eqs.([l5|) and ( |l6| ) and compare our simula- 
tions with the theoretical prediction from eqs.(|l2|) and 
(|l3|). The quantitative coincidence is within a few per- 
cent. This result is even more convincing when keeping 
in mind that the theoretical prediction is just a first or- 
der computation in the control amplitude which actually 
is not so small. For larger values of K we observe an ad- 
ditional frequency in the control signal. Such a property 



can be attributed to the second (stable) Floquet branch 
and may also be evaluated by eq.(0). Finally, as a by- 
product, we obtain an estimate of the Floquet exponent 
A + iuj of the uncontrolled orbit. 




2n 



0.7 
0.5 
0.3 







i i ^^^A^yyy^ 




000000 





2.5 



K 



3.5 



FIG. 2. Real and imaginary part of the Floquet exponent 
for a numerical simulation of eqs.(^). The data (symbols) 
have been obtained from the decay of the control signal on 
variation of the control amplitude K. The solid lines indicate 
the analytical solutions of eqs. ([if) and © with A w 0.97, 
X ' « -0.23. 

In addition to analytical calculations and computer 
simulations we have performed experiments on a non- 
linear electronic circuit (see e.g. 10 1). We consider a 
nonlinear diode resonator consisting of a capacity diode 
(1N4005), an inductor (470/^H), and a resistor (40fi) 
(cf. fig.|. 




FIG. 3. Experimental setup of the nonlinear diode res- 
onator with time-delayed feedback device. 

The control device consists of a cascade of electronic 
delay lines with a limiting frequency of about 2MHz and 
several operational amplifiers acting as preamplifier, sub- 
tractor, or inverter. The device allows to apply a con- 
trolling force of the form ±K[U(t) - eU(t - r)] + U 
with parameter ranges K = . . . 300, e — ... 2, r = 
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10 . . . 7000ns, and Uq = —5 . . . + 5V. For conventional 
delayed feedback control e has carefully to be adjusted to 
one and the offset Uq to zero. This was done in the exper- 
iment reported here. The circuit was sinusoidally driven 
with an amplitude of 2V and a frequency of 990kHz. Ac- 
cordingly the delay time was set to r = 1010ns. From the 
transient dynamics of the control signal we again obtain 
the decay rate and the frequency (cf. fig-^)- For compar- 
ison with the theoretical prediction we take the values 
K c ks 34 and K opt ~ 37.5 to determine the two parame- 
ters \'t and At from eqs.( |T5"|) a nd (|l6|). The quantitative 
agreement with eqs.(p"2[), (|13|) is within a few percent, 
except for the real part beyond K opt . 
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FIG. 4. Real and imaginary part of the Floquet exponent 
for the electronic circuit experiment. The data (symbols) have 
been obtained from the decay of the control signal on vari- 
ation of the control amplitude K. The solid lines indicate 
the analytical solutions of eqs,(|l2|) and ( |l3| ) with Ar ~ 1.09, 
X'r w -0.016. 

Apart from the reasons already mentioned the devi- 
ations can be attributed to the limited accuracy of the 
value K opt . Since the transients are affected by noise a 
precise estimate of the exponents is difficult to obtain for 
small decay rate A. 

IV. CONCLUSION 

We have shown that the main limiting factor for time 
delayed feedback control results from the torsion of the 
unstable periodic orbit. This topological property de- 
termines whether the control mechanism works at all. 
We have worked out the general features of the transient 
behaviour including critical and optimal control ampli- 
tudes. Our approach describes at least the generic prop- 
erties for stabilizing unstable periodic orbits with an un- 
stable manifold like a Mobius strip. Our simulations 
have shown that the features described above are acces- 



sible from the transient behaviour of the control signal 
and hence are observable in experiments. The electronic 
circuit experiment demonstrates that an analysis along 
these lines is possible even for ultrafast experiments. Our 
theoretical approach, based on an expansion of the gen- 
eral expression (J?]), resembles a Ginzburg-Landau like 
treatment of phase transitions. It can be easily extended 
to incorporate e.g. the degeneracy of several Floquet ex- 
ponents or the features of spatially extended, that means 
high-dimensional systems. 
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